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Abstract 

Generating high quality, structured, continuous, 
body-fitted grid systems (multiblock grid systems) for 
complicated geometries has long been a most labor- 
intensive and frustrating part of simulating flows in 
complicated geometries. Recently, new methodologies 
and software have emerged that greatly reduce the hu- 
man effort required to generate high quality multiblock 
grid systems for complicated geometries. These meth- 
ods and software require minimal input from the user — 
typically, only information about the topology of the 
block structure and number of grid points. This paper 
demonstrates the use of the new breed of multiblock 
grid systems in simulations of internal flows in compli- 
cated geometries. The geometry used in this study is 
a duct with a sudden expansion, a partition and an ar- 
ray of cylindrical pins. This geometry has many of the 
features typical of internal coolant passages in turbine 
blades. The grid system used in this study was gen- 
erated using a commercially available grid generator. 
The simulations were done using a recently developed 
flow solver, TRAF3D.MB, that was specially designed 
to use multiblock grid systems. 

Introduction 

Generating high quality, structured, continuous, 
body-fitted grid systems (multiblock grid systems) for 
complicated geometries has long been a most labor- 
intensive and frustrating part of simulating flows in 
complicated geometries. Although the grid generation 
techniques that have been used are often termed “auto- 
matic,” 1 they have in fact been semi-automatic at best 
since considerable human involvement has been required 
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to generate acceptable grids. The human involvement 
has been to “carve up” the flow domain into zones 
of reasonable shape and, through trial and error, set 
various parameters for the grid generators that opti- 
mize the grid system with respect to the conflicting 
objectives of (a) clustering grid points where they are 
needed, (b) maximizing “smoothness” (e.g., minimize 
rapid changes in grid spacing), and (c) minimizing skew- 
ness. To avoid the difficulty in generating body-fitted, 
structured grid systems, researchers have invented dif- 
ferent approaches to the “gridding” of flow domains. 
Using separate body-fitted grids for different parts of 
the geometries and letting the component grids overlap 
in an arbitrary manner (Chimera grids 2-4 ), has been a 
popular choice (e.g., Ref’s 5-7). The use of unstruc- 
tured grids has also received much attention in recent 
years, precisely for the relative automation of the grid 
generation process for complex geometries and for the 
apparent flexibility that unstructured grids offer to dy- 
namically adapt the grid to the computed solutions. 

Recently, new methodologies and new software 8-12 
have emerged that greatly reduce the human effort re- 
quired to generate high quality multiblock grid systems 
for complicated geometries. These methods and soft- 
ware require minimal input from the user — typically, 
only information about the topology of the block struc- 
ture and density of grid points. The grid generation 
software then optimizes the grid with respect to orthog- 
onality, smoothness and other conflicting objectives. A 
distinguishing feature of the multiblock grid systems 
is that they, in general, cannot be mapped to a sin- 
gle index (or computational) space. In essence, these 
multiblock grid systems are unstructured at the block 
level but structured at the cell level (array arrangement 
of cells within each block). The absence of stucture 
at the block level results in great flexibility in the cre- 
ation of mesh topologies for complicated geometries and 
an opportunity to generate high quality grid systems. 
The only requirement for the block structure is that the 
blocks be topologically quadrilateral in two dimensions 
(2D) and hexahedra in three dimensions (3D). 
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The objective of this paper is to demonstrate the 
use of the new breed of multiblock grid systems in simu- 
lations of internal flows in complicated geometries. The 
geometry used in this study is a duct with a sudden ex- 
pansion, a partition and an array of cylindrical pins 
(see Fig. 1). This geometry has many of the features 
of internal coolant passages in turbine blades and has 
been used in experiments aimed at increasing the un- 
derstanding of coolant flow and heat transfer in turbine 
blades 1 **. The grid system used in this study was gen- 
erated using a commercially available grid generator 
called GridPro. 12 The simulations were done using a 
recently developed flow solver, TRAF3D.MB, that was 
specially designed to use multiblock grid systems 14 . 

The paper is organized as follows: Following this 
introduction, the problem to be studied is described. 
Afterwards, the grid system used in the computations 
is discussed and some important features pointed out. 
Next, the major elements of flow' solver will be briefly 
described. Finally, results of computations are pre- 
sented. The paper ends with some concluding remarks. 

Description of Problem 

The problem to be analyzed is the flow of air 
through a duct with a sudden expansion and, down- 
stream of the expansion, a partition that splits the duct 
in two branches and an array of cylindrical pins. A 
schematic drawing of the duct is shown in Fig. 1. As 
indicated in the figure, the air flows through the duct 
from left to right. The array of pins below the partition 
creates a blockage to the flow in the lower branch of the 
duct that acts to divert some of the flow through the 
upper branch. The floor of the duct is heated. 

The flow field in the branched duct is quite com- 
plicated. Nonetheless, some qualitative features of the 
flow are predictable based on experience and intuition. 
In particular, separation of the flow at the expansion 
corner and on the upper side of the partition is ex- 
pected. However, the location of reattachment points 
on the upper side wall and partition, the extent of the 
recirculation region above the partition, and the effect 
of the pin array on the flow field cannot be predicted 
easily without numerical simulations. Consequently, 
heat transfer and wall temperature patterns in the duct 
cannot be reliably predicted without detailed numerical 
investigations. 

The branched duct geometry used in the current 
study has many of the features found in coolant pas- 
sages of turbine blades, such as the pins and the sud- 
den expansion. Also, the heating of the floor emulates 
the thermal loading on a turbine blade. The duct was 
designed and used in experiments at NASA Lewis Re- 


search Center 13 in a project aimed at enhancing un- 
derstanding of flows in turbine blade coolant passages. 

In the experiments, measurements of flow field and 
heat transfer were done for inflow Reynolds numbers 
of 45,000, 335,000, and 726,000 (based on flow speed 
in the entrance section and hydraulic diameter of the 
inflow duct), corresponding to Mach number from 0.03 
to 0.45. Heat transfer results from the experiments 
were reported in Ref. 13. 

Grid System 

A multiblock grid system was generated for the 
branched duct geometry shown in Fig. 1. The grid sys- 
tem was generated using a commercially available grid 
generator called GridPro. 12 A 2-D grid system was ini- 
tially generated for the cross section shown in Fig. lb. 

A stack of the 2-D grids forms the final 3-D grid sys- 
tem, which extends from the channel floor to the sym- 
metry plane midway between the floor and the ceiling. 
Figures 2-4 show the 2-D grid system used in the cur- 
rent study. Figure 2a shows the block structure of the 
grid whereas Fig. 2b shows a close-up view of the block 
structure around the partition and pins. Figure 3 shows 
grid system for the main-part of the duct and Fig. 4 
shows a close-up view of the grid near the corner at the 
sudden expansion and around the partition and pins. 

There are several features to note about the grid 
system used in this study. First, as seen in Fig. 2a, 
there are several singular points in the grid where three, 
five or six blocks come together. Because of these sin- 
gular points, the grid system cannot be mapped to a 
single computational space, and flow solvers utilizing 
grid systems of this type must be set up to handle the 
singularities properly. However, allowing these kinds of 
singularities in the grid is exactly what gives the flexi- 
bility that is needed for generating high quality contin- 
uous structured grids for complex geometries such as 
the branched duct. A second feature to note is that 
high quality “boundary-layer type” grids are obtained 
around all solid surfaces. Thus, O-grids are formed 
around the pins and grid lines wrap around the par- 
tition (see Fig. 4a). Also, grid fines run parallel to 
the side walls of the duct and turn around the slightly 
rounded expansion corner (note, in the actual geom- 
etry used in the experiments reported in Ref. 13, the 
expansion corner is slightly rounded, here modeled with 
a radius of 0.5mm). These high quality boundary-layer 
type grids lend themselves well to computations of vis- 
cous flows and facilitates accurate implementation of 
physical boundary conditions. A final feature to note 
about the grid system in Fig. 2-4 is that the grid every- 
where smooth and, except in the immediate vicinity of 
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singularities, nearly orthogonal. In particular, the grid 
is very nearly orthogonal at all solid surfaces, as the 
singularities are located away from all physical bound- 
aries. 

The 3-D grid system used in the current study 
contains about 2 million cells. Each 2-D cross section 
of the grid contains about 50,000 grid points, 65% of 
which are in the blocks immediately adjacent to solid 
wall surfaces (the boundary layer grids). The grid sys- 
tem was generated in 283 full-face matching blocks. 
The size of the 2-D blocks varies from 8 X 8 to 65 x 24. 
Typical grid spacing normal to sidewalls, pins and par- 
tition is 2 • 10 -4 whereas grid spacing at the channel 
floor is 1.1 * 10”* (note, dimensions are normalized by 
the height from the floor of the channel to the symme- 
try plane). 

Flow Solver 

The simulations performed in this study were done 
using a recently developed computer code called TRAF- 
3D.MB. 14 This code is is a general purpose flow solver, 
designed specifically for use in simulations of flows in 
complicated geometries. The code is based on an earlier 
code, TRAF3D, an efficient code designed for simula- 
tions of flows in turbine cascades. 15 The TRAF3D.MB 
code uses a cell-ceneterd finite volume discretization 
of the compressible Navier-Stokes equations. The dis- 
cretization is second-order accurate in space. The code 
offeres a choice between central differencing with artifi- 
cial dissipation and the ATJSM 16 scheme for the inviscid 
flux terms of the governing equations. Central differ- 
encing is used for diffusion terms. Baldwin-Lomax eddy 
viscosity model 17 ' 15 is currently used to model the ef- 
fects of turbulence on the mean flow. The model is 
applied at all solid walls of the geometry. In corners 
between two walls, the contributions to the eddy vis- 
cosity from the walls is averaged as shown in Ref. 15. 
Boundary conditions are implemented using ghost cells. 
The types of boundary condition to apply at physical 
boundaries are conveniently specified at runtime. 

Solutions are marched to steady state using four 
stage Runge-Kutta time stepping scheme. Convergence 
to steady state is accelerated by using local time step- 
ping, implicit residual smoothing and multigrid. To re- 
duce overhead due to communication of data between 
blocks in multiblock grid systems, exchange of data 
between blocks takes place only before the first stage 
of the Runge-Kutta scheme. Also, implicit residual 
smoothing is done independently within each block. 

The TRAF3D.MB code was specifically designed 
for use with mutliblock grid systems such as the one 
used in the current study. The code is eqquipped with 


a completely general multiblock capability, i.e., there 
is no limitation on the number of blocks a grid system 
may have or on how they can be connected. Ghost 
cells are used to accomplish the necessary communica- 
tion between the blocks of the multiblock grid system. 
The ghost cells, in effect, create an overlap between 
blocks and communication between blocks simply in- 
volves loading into the appropriate ghost cells the data, 
needed from neighboring blocks. In the current imple- 
mentation two layers of ghost cells are used around each 
block of the grid system, creating, in effect, a four-cell 
wide overlap between blocks. This overlap allows the 
use of fourth-difference artificial dissipation or up to 
third order accurate upwind formulas for the inviscid 
terms of the governing equations without loss of accu- 
racy at blocks boundaries. Communication takes place 
on all levels of the multigrid cycle at the beginning of 
each smoothing sweep. 

The TRAF3D.MB code has been used in studies 
of various turbomachinery and prototype flows. Two 
simulations in particular serve as validations of the ac- 
curacy and capability of the flow solver. The first was 
a simulation of a laminar flow over a backward facing 
step. Exellent agreement with available experimental 
data was obtained for both reattachement length and 
velocity profiles, with reattachment length predicted to 
within 3% of the experimentally measured values. The 
results are reported in Ref. 18. The second simulation 
was of the flow and heat transfer in a turbine rotor 
cascade, where detailed resolution of the flow field in 
the tip clearance region was obtained. The predicted 
values of heat transfer on the rotor tip matched very 
well with experimentally determined values. The re- 
sults frome these numerical simulations are described 
in Ref. 19. 

Results 

Simulations of flow in the branched duct are cur- 
rently under way. At the time of printing of this paper, 
only very preliminary results are available. The results 
presented here are for flow conditions modeled after the 
medium Reynolds number case (Re = 335000) reported 
in Ref. 13. The appropriate flow conditions were ob- 
tained by specifying, along with the Reynolds number, 
the static back pressure at the exit of the lower and 
upper branches of the duct to be the same as in the 
experiments (pressure specified as a ratio of the inlet 
total pressure). 

The boundary conditions used for the computa- 
tions presented here are as follows: The static pressure 
ratios at the exit of the upper and lower branches of 
the duct were specified as 0.9685 and 0.9725, respec- 
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tively. At the inlet, a constant total pressure was spec- 
ified, equivalent to assuming a flat incoming velocity 
profile. At all solid walls, a no-slip condition was im- 
posed. Wall temperature was fixed at 1.1 times the inlet 
total temperature. These boundary conditions are only 
an approximation to the actual run conditions used in 
the experiments. The most significant difference is that 
in the experiments, a flow with a developed boundary 
layer entered the test section. Consequently, the mass 
flow rate obtained in the simulations is likely in excess 
of the experimental value. 

A sample of results from the simulations is shown 
in Fig’s 5-8. Figure 5 shows the velocity field in the 
symmetry plane of the duct. Only shown is the region 
around the pin array and partition, where the topol- 
ogy of the grid is most complicated. Note, velocity 
at every-other grid point is plotted. Figure 6 shows a 
close up view of the flow around some of the pins and 
the leading edge of the partition. Figure 5 reveals how 
the array of pins acts as a blockage that diverts the flow 
towards the upper channel of the duct. Consequently 
the flow approaches the partition at an angle of attack 
close to 60°. The flow separates above the partition 
as expected. The blockage effect also causes the flow 
to enter the pin array with a non-uniform angle. This 
contributes to a complicated flow field within the pin 
array (see e.g., Fig. 6) and behind the pins. 

Figure 7 shows the velocity field near the expan- 
sion corner. As in Fig. 5 and 6, the flow in the sym- 
metry plane of the duct is plotted. Every other vector 
on grid lines normal to the solid wall is plotted. As the 
figure shows, the flow separates at the corner. The de- 
tails of the flow at the corner are well captured without 
excessive resolution away from the corner. 

Finally, Fig. 8 shows contours of static pressure 
ratio around the leading edge of the partition. Note 
the continuity of the contours across all block bound- 
aries. As a rule, this continuity is obtained when multi- 
block grids are used since no interpolations are needed 
at block boundaries. 

Conclusions 

Although only preliminary computations of the 
branched duct flow have been formed thus far, it is 
still possible to draw important conclusions about the 
suitability of the present methodology for simulation of 
viscous flows in complex geometries. A first conclusion 
is related to the use of multiblock grid systems to model 
the complicated geometry. A high quality multiblock 
grid system was successfully generated for the compli- 
cated branched duct geometry. The grid obtained was 
smooth and nearly orthogonal. In particular, high qual- 


ity grid was obtained around all solid wails. In the 2-D 
cross section, about 65% of grid points were located in 
the blocks around the pins, partition and side walls, 
with the remaining 35% located away from the walls. 
Thus, one may conclude that other types of structured 
grids, such as Chimera grids, with similar resolution of 
the boundary regions, would require about the same 
total number of grid points as the current grid scheme. 
On the other hand, the use of complicated interpo- 
lation schemes to transfer data between blocks, as is 
needed when Chimera grids are used, is avoided in the 
the multiblock flow solver. 

In the current approach, the grid system was gen- 
erated in many small blocks. Although it is possible to 
use fewer blocks to generate the grid, the fact remains 
that multiblock grid systems for complicated geome- 
tries tend to come in many small blocks. For two rea- 
sons, this can be a drawback. First, for flow solvers that 
run on vector supercomputers, the small blocks tend to 
lead to short vector lengths and degradation in per- 
formance. This effect was observed for the flow solver 
used in this study. This can be remedied by merging 
small blocks and form larger ones that allow greater 
vector lengths, or by writing the computer codes in 
special ways that allow greater vector lengths. Second, 
using many small blocks will impact negatively the rate 
of convergence to steady state if coupling between the 
blocks is not sufficiently strong. At this time definitive 
conclusions about the importance of this effect as it re- 
lates to the present flowsolver cannot be drawn as the 
needed numerical experiments have not yet been per- 
formed. Indications are, however, this is not a serious 
effect and that coupling between blocks in the present 
algorithm is strong enough. 

A final observation to be made relates to the use 
of Baldwin-Lomax turbulence model. Although the 
Baldwin-Lomax model has been found to work surpris- 
ingly well in many complicated flows (see, e.g., Ref. 19) 
it can not be expected to work perfectly in the current 
application, due to the complicated 3-D flow field, for 
instance behind the array of pins in the lower branch 
of the duct. For future detailed simulations, a more 
general two-equation turbulence model should be used. 
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Note: Passage height Is 5.08 
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Figure 1. Branched duct geometry (a) experimental setup (b) schematic drawing of the test section 
(all dimensions are in centimeters). Figure reprinted from Ref. 13 (with permission). 
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Figure 2. Topology of multiblock grid system (a) overall grid, (b) around pin array and par- 
tition 



Figure 3. Grid system for the center region of the branched duct. 
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(a) 



(b) 

Figure 4. Close-up view of grid near (a) expansion comer 
and (b) pin array and position 
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Figure 5. Velocity field in the lower branch (symmetry plane) 



(O ») 

Figure 6. Close-up view of velocity field near (a) pin 3, (b) pin 7, (c) 

pin 12 and (d) leading edge of the partition. 


9 













Form Approved 

REPORT DOCUMENTATION PAGE ombno. 0704 -om 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of Information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 

1. AGENCY USE ONLY ( Leave blank) I 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 


November 1995 Contractor Report 


4. TITLE AND SUBTITLE 

Computations of Viscous Flows in Complex Geometries Using Multiblock 
Grid Systems 

5. FUNDING NUMBERS 

WU-505-90-5K 

6. AUTHOR(S) 

Erlendur Steinthorsson and Ali A. Ameri 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

Institute for Computational Mechanics in Propulsion 
22800 Cedar Point Road 
Cleveland, Ohio 44142 

E-9992 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

National Aeronautics and Space Administration 
Washington, D.C. 20546-0001 

. . i 

NASA CR-198419 
ICOMP 95-23 
AIAA-95-0177 

11. SUPPLEMENTARY NOTES 

Prepared for the 33rd Aerospace Sciences Meeting and Exhibit sponsored by the American Institute of Aeronautics and Astronautics, Reno, Nevada, 
January 9-12, 1995. Erlendur Steinthorsson, Institute for Computational Mechanics in Propulsion, NASA Lewis Research Center, work funded under 
NASA Cooperative Agreement NCC3-370); Ali A. Ameri, AYT Corporation, Beachwood, Ohio. ICOMP Program Director, Louis A. Povinelli, 
organization code 2600. (216) 433-5818. 

12a. Dl STRIBUTIO N/A V Al L AB ILIT Y STATEMENT 

12b. DISTRIBUTION CODE 

Unclassified - Unlimited 
Subject Categories 34 and 64 


This publication is available from the NASA Center for Aerospace Information, (301) 621-0390. 


13. ABSTRACT (Maximum 200 words) 


Generating high quality, structured, continuous, body-fitted grid systems (multiblock grid systems) for complicated 
geometries has long been a most labor-intensive and frustrating part of simulating flows in complicated geometries. 
Recently, new methodologies and software have emerged that greatly reduce the human effort required to generate high 
quality multiblock grid systems for complicated geometries. These methods and software require minimal input from the 
user-typically, only information about the topology of the block structure and number of grid points. This paper demon- 
strates the use of the new breed of multiblock grid systems in simulations of internal flows in complicated geometries. The 
geometry used in this study is a duct with a sudden expansion, a partition and an array of cylindrical pins. This geometry 
has many of the features typical of internal coolant passages in turbine blades. The grid system used in this study was 
generated using a commercially available grid generator. The simulations were done using a recently developed flow 
solver, TRAF3D.MB, that was specially designed to use multiblock grid systems. 


14. SUBJECT TERMS 

CFD; Multigrid; Multiblock grid; Turbine cooling 

15. NUMBER OF PAGES 
12 

16. PRICE CODE 

A03 

17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 



18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 

19. SECURITY CLASSIFICATION 
OF ABSTRACT 

Unclassified 

20. LIMITATION OF ABSTRACT 


NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 


Prescribed by ANSI Std. Z39-18 
298-102 

















